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ABSTRACT 

A  coupled  atmosphere-ocean  model  developed  at  the  Institute  for  Space  Studies  at  NASA  Goddard 
Space  Flight  Center  (Russell  et  al.,  1995)  was  used  to  verify  the  validity  of  Haney-type  surface  thermal 
boundary  condition,  which  linearly  connects  net  downward  surface  heat  flux  Q  to  air  /  sea  temperature  dif¬ 
ference  NT  by  a  relaxation  coefficient  k.  The  model  was  initiated  from  the  National  Centers  for 
Environmental  Prediction  (NCEP)  atmospheric  observations  for  1  December  1977,  and  from  the  National 
Ocean  Data  Center  (NODC)  global  climatological  mean  December  temperature  and  salinity  fields  at  1°  x 
1°  resolution.  The  time  step  is  7.5  minutes.  We  integrated  the  model  for  450  days  and  obtained  a  complete 
model-generated  global  data  set  of  daily  mean  downward  net  surface  flux  Q,  surface  air  temperature  T A , 
and  sea  surface  temperature  T0.  Then,  we  calculated  the  cross-correlation  coefficients  (CCC)  between 
Q  and  AT,  The  ensemble  mean  CCC  fields  show  (a)  no  correlation  between  Q  and  AT  in  the  equatorial  re¬ 
gions,  and  (b)  evident  correlation  (CCC^  0.7)  between  Q  and  AT  in  the  middle  and  high  latitudes. 
Additionally,  we  did  the  variance  analysis  and  found  that  when  k—  120  W  m-2K-1,  the  two  standard  devia¬ 
tions,  Oq  and  <7KAr ,  are  quite  close  in  the  middle  and  high  latitudes.  These  results  agree  quite  well  with  a 
previous  research  (Chu  et  aL,  1998)  on  analyzing  the  NCEP  re-analyzed  surface  data,  except  that  a  smaller 
value  of  k  (80  W  m_2K-1)  was  found  in  the  previous  study. 

Key  word®  Air— sea  coupled  system.  Ocean  surface  fluxes,  Surface  thermal  boundary  condition 

1.  Introduction 

The  ocean  model  is  thermally  driven  by  the  net  downward  flux  of  heat  across  the  ocean 
surface.  Q,  which  is  the  sum  of  the  downward  flux  of  solar  radiation,  Rs ,  minus  the  net  up¬ 
ward  flux  of  longwave  (or  “back”  radiation),  Rb  ,  sensible  heat  Q  H ,  and  latent  heat  Q  E , 

Q=  Rs-(Rb+Q„+QE)  •  U) 

Under  the  assumption  that  the  ocean  is  in  contact  with  an  atmospheric  equilibrium  state, 
Haney  (1971)  obtained  a  very  simple  heat  flux  formulation, 

2=K(r;-r0)  (2) 

where  T'k  is  the  apparent  atmospheric  equilibrium  temperature,  T0  the  sea  surface 
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temperature  (SST),  and  k  is  a  relaxation  coefficient.  In  spite  of  its  temporal  variation,  the 
parameter  is  usually  taken  as  positive  values  between  10  to  50  W  m-2K-1  by  various  authors 
(e.g.,  Marotzke  and  Stone,  1995;  Cai  and  Chu,  1996).  Here  Tk  should  be  computed  from 
surface  heat  fluxes  and  their  dependence  on  temperature,  as  done  by  Han  (1984)  and 
Oberhuber  (1988).  Despite  Haney's  emphasis  on  that  Q  is  proportional  t  ofj  -  T0,  many 
subsequent  investigators  have  replaced  T'k  by  T  A .  Here,  T  A  is  the  surface  air  temperature 
(SAT).  Thus,  we  call 

2  =  kAT,  AT  =  T  A  —  T0  (3) 

the  Haney-type  surface  boundary  condition.  Another  trend  is  to  separate  thesolar  radiation 
from  the  net  downward  flux  (e.g.,  Blumberg  and  Mellor,  1987) 

Q=RS+Q,  Q=-(Rb+QH  +  Q*)  (4) 

and  to  replace  Q  by  Q  in  the  Haney-type  surface  boundary  condition 

2=  kAT  .  (5) 

We  call  2  the  net  downward  heat  flux  with  solar  radiation,  and  Q  the  downward  heat  flux 
without  solar  radiation. 

Chu  et  al.  (1998)  used  the  global  reanalyzed  data  of  Q,  surface  air  temperature  T A ,  and 
sea  surface  temperature  T0  from  the  National  Centers  for  Environmental  Prediction  (NCEP) 
from  1  October  1994  to  31  December  1995  to  verify  the  validity  of  Haney-type  surface 
thermal  boundary  conditions.  They  found  that  the  Haney-type  surface  thermal  condition  (2) 
represents  the  net  heat  flux  quite  well  for  the  middle  and  high  latitudes,  but  very  poor  in  the 
equatorial  regions.  Is  the  result  universal  or  model— dependent?  We  should  use  a  different  type 
model  to  verify.  Since  NCEP  constantly  has  available  observational  data  assimilation,  we 
chose  a  coupled  air— ocean  model  without  data  assimilation  in  this  study  to  see  if  the  previous 
result  (Chu  et  al„  1998)  still  holds. 

2.  A  coupled  atmospheric-ocean  model 

2. 1  M  odel  description 

A  coupled  atmosphere— ocean  model  (Russell  et  al.,  1995)  developed  at  the  NASA 
Goddard  Institute  for  Space  Studies  (GISS)  was  used  for  this  study.  The  atmospheric  model  is 
similar  to  that  of  Hansen  et  al.  (1983)  except  that  the  atmospheric  dynamic  equations  for 
mass  and  momentum  are  solved  using  C  grid  scheme  and  the  advection  of  potential  enthalpy 
and  water  vapor  uses  the  linear  upstream  scheme  (Russell  and  Lerner,  1981).  The  global 
ocean  model  conserves  mass,  allows  for  divergent  flow,  has  a  free  surface  and  uses  the  linear 
upstream  scheme  for  the  advection  of  potential  enthalpy  and  salt,  Both  models  run  at  4°  x 
5°  resolution,  with  9  vertical  layers  for  the  atmosphere  and  13  layers  for  the  ocean.  Twelve 
straits  are  included,  allowing  for  subgrid-scale  water  flow.  Runoff  from  land  is  routed  into 
approximate  ocean  basins.  Atmospheric  and  oceanic  surface  fluxes  of  water,  heat  (excluding 
solar  radiation),  and  momentum  are  of  opposite  sign  and  are  applied  synchronously.  Flux  ad¬ 
justments  are  not  used.  Except  for  partial  strength  alternating  binomial  filters  (Shapiro, 
1970),  which  are  applied  to  the  momentum  components  in  the  atmosphere  and  oceans,  there 
is  no  explicit  horizontal  diffusion.  For  more  information  about  this  coupled  model,  readers 
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are  referred  10  Russell  el  ul  11995).  The  GISS  coupled  model  was  initialed  from  the  NCEP 
atmospheric  observations  lor  I  December  1977.  and  from  the  National  Ocean  Data  Center 
(NODC)  global  climatological  mean  December  temperature  and  salinity  fields  with  l°* 
1°  resolution  The  time  step  is  7,5  minutes,  and  the  model  was  integrated  for  450  days.  Our 
purpose  here  is  to  build  up  a  dynamically  and  thermodynamically  consistent  data  set  of 
global  surface  fluxes  Rs  ,/ih.£n  ,Qr ,  and  temperatures  T0,T  \  .  rather  than  to  study  physical 
processes  Thus,  a  450-day  time  integration  will  provide  sufficient  data  to  fulfill  this  task 
However,  we  should  bear  in  mind  that  the  model  surface  conditions  must  change  from  year  to 
year  so  that  the  results  in  this  study  may  be  modified  quantitatively, 


Longitude  (E).  Depth=0(m) 
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Pi*  I  Monthly  mean  sea  surface  temper, nine  <"C)  simulated  bythc  CUSS  Coupled  Model.  U1 
January  I97H.  «nd  lb)  July  |97ji 
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2.2  Model  output 

We  choose  model  generated  global  surface  fluxes  Rs  ,Rb ,QH,QE,  and  temperatures 
rG,rA  data  on  4°x  5°  grid  for  this  study.  Since  Haney-type  surface  thermal  boundary  con¬ 
ditions  are  commonly  used  for  ocean  modeling  on  a  seasonal  time  scale  or  longer,  we  first 
computed  daily  means  of  these  variables  to  filter  out  high  frequency  variability  and  to  obtain 
a  data  set  of  global  Q(x,y,d),Q(x,y,d),TA(x,y,d),T0  (xj,d)  for  each  day,  d,  from  1  February 
1978  to  5  May  1979. 

2.2.1  Sea  surface  temperature 

The  GISS  coupled  model  simulates  a  number  of  interesting  sea  surface  temperature 
(SST)  patterns,  indicative  of  the  interaction  between  ocean  and  atmosphere  (Fig.  1).  Western 
boundary  currents  such  as  the  Gulf  Stream,  Kuroshio,  Brazil,  and  Agulhas  currents  off  the 
eastern  shores  of  North  America,  Japan,  South  America,  and  Southern  Africa,  respectively 
appears  as  meandering  bands  of  warm  waters  moving  poleward  into  the  ocean  interior.  Large 
scale,  seasonal  variation  is  also  simulated.  During  January  (Fig.  la)  the  Southern  Hemisphere 
is  warm  while  the  cooler  temperatures  are  located  in  the  northern  temperate  zone.  The  con¬ 
verse  is  true  in  July  (Fig.  lb)  where  we  see  the  other  phase  of  the  seasonal  cycle.  Upwelling 
and  advection  leads  to  large  asymmetries  in  the  subtropical  thermal  structure.  The  simulated 
SST  fields  show  the  complex  connections  of  the  California  Current  with  the  equatorial  cur¬ 
rent  system.  South  Pacific  eastern  boundary  current  and  Peru  current. 

Furthermore,  the  model  simulates  the  asymmetric  distribution  of  warm-cool  waters  in 
the  equatorial  regions.  We  see  that  the  warm  water  is  not  only  within  ocean  basins,  but  also 
globally  with  the  large  warm  pool  in  the  western  Pacific,  Indonesia,  and  eastern  Indian 
Ocean.  The  cool  water  is  located  in  the  eastern  Pacific.  This  leads  to  asymmetric  along-equa- 
torial  atmospheric  circulation  -  the  Walker  Cells. 

2.2.2  Surface  air  temperature 

The  GISS  coupled  model  simulates  realistic  surface  air  temperature  (SAT)  variabilities. 
On  monthly  fields,  a  general  parallelism  exists  between  SAT  and  SST  (Fig.  2).  Large  scale, 
seasonal  variation  in  SAT  is  also  simulated.  During  January  (Fig.  2a)  the  Southern  Hemi¬ 
sphere  is  warm  while  the  cooler  temperatures  are  located  in  the  northern  temperate  zone.  The 
converse  is  true  in  July  (Fig.  2b)  where  we  see  the  other  phase  of  the  seasonal  cycle.  Large 
scale  zonal  asymmetries  were  also  simulated.  SAT  is  low  over  the  upwelling  regions,  such  as 
California  coast.  Furthermore,  the  model  simulates  the  asymmetric  distribution  of 
warm— cool  SAT  in  the  equatorial  regions.  We  see  that  the  SAT  is  warmer  over  the  western 
Pacific  warm  pool  and  SAT  is  cool  over  the  eastern  Pacific. 

2.2.3  Ocean-atmosphere  heat  exchange 

The  GISS  coupled  model  simulates  realistic  variabilities  of  heat  budget  at  the  air— ocean 
interface,  including  radiation,  sensible  heat  flux,  and  latent  heat  flux. 

Short-wave  radiation 

The  amount  of  solar  radiation  actually  incident  on  the  ocean  surface  depends  on  the  time 
of  year,  the  time  of  day,  the  latitude,  and  the  atmospheric  absorption  and  scattering  (especial¬ 
ly  by  clouds.)  Solar  radiation  is  either  absorbed  by  the  sea  or  reflected  from  its  surface.  The 
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Fig  2  Monihly  mean  surface  air  temperature  t"C)  simulated  bythe  OISS  Coupled  Model:  U) 
January  1978,  and  (h)  July  1978, 


degree  of  reflection  (or  called  Ihe  ocean  albedo)  used  by  the  GISS  coupled  model  is  primary  a 
function  of  solar  zenith  angle  (.0)  and  surface  wind  speed  (F, ): 

a.  =  0.012+  0.042  l.r:  +  0.1283*J  -  0.04.V4  +  ;  . 

where  v“  I  cos  '0,  This  parameterization  is  based  on  calculations  of  Fresnel  reflection 
from  Ihe  ocean  surface  as  a  function  of  V t  as  specified  by  Cox  and  Munk  (1956).  Figure  3 
shows  the  GISS  model  simulated  seasonal  variation  of  the  net  downward  solar  radiation  over 
the  global  ocean  surface  (fls  ),  In  January  (July),  the  Southern  (Northern)  Hemisphere  ocean 
surface  receives  more  solar  energy  than  the  other  hemisphere.  For  example,  an  area  enclosed 
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SW  flux  monthly  mean  Jan  1978 


SW  flux  monthly  mean  Jul  1978 


Fig  J  Monthly  mean  surface  tie l  downward  solar  raduitinn  (W  nr*  simulated  by  the  CUSS 
Coupled  Model  In)  January  IU7K.  and  ihl  July  IU7& 

b>  the  250  W  nr  contour  in  the  central  Pacific,  shifted  toward  the  Southern  (Northern  1 
Hemisphere  in  January  (July)  1978, 

Long-wave  radiation 

Besides  absorbing  and  reflecting  solar  radiation,  the  ocean  surface  also  emits  radiation  ol 
a  wavelength  appropriate  to  its  temperature  and  absorbs  long-wave  radiation  transmitted 
downward  from  clouds  and  the  atmosphere.  Figures  4a  and  4b  show  the  G1SS  coupled  model 
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LW  flux  monthly  mean  Jan  1978 


Longitude  (E),  Depth=0(m) 


LW  flux  monthly  mean  Jul  1978 


Longitude  (E),  Depth=0(m) 


Fi(.  A  Monthly  mcun  surlucc  net  upward  long  wave  railiaiion  i W  m’>  simulated  by  Ihc  C.ISS 
Coupled  Model;  (a)  January  1978,  and  (bt  July  1 M7K 


outpul  of  mean  January  and  July  <  19~8 )  ncl  upward  long-wave  radiative  Dux  at  the  ocean 
surface  </?b  1.  respectively  The  magnitude  of  is  much  smaller  than  the  magnitude  of  Ks . 
which  makes  the  radiation  balance  of  the  ocean  surface  is  everywhere  positive.  This  agrees 
quite  well  with  earlier  sludics  (e.g..  Budyko.  19561, 

Sensible  heat  flux 


Figures  5a  and  5b  show'  the  GISS  coupled  model  output  of  mean  January  and  July 
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SH  (lux  monthly  mean  Jan  1978 


Longitude  (E), 
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SH  flux  monthly  mean  Jul  1978 


Longitude  (E), 


-30 


-60 

30 


hg  5  Monthly  mean  surface  sensible  heat  tlnx  tW  nr '  umulurcd  bythc  CUSS  Coupled  Model 
u  I  January  1978,  and  (b)  July  1978, 


(1978)  upward  sensible  heat  flux  at  the  ocean  surface  (£H  ).  respectively,  Large  values  of 
cun  be  expected  in  air  which  is  very  much  cooler  lhan  the  surface  of  the  ocean  across  w  hich  it 
is  flowing,  especially  when  the  wind  is  strong  Therefore,  il  L%  not  surprising  to  find  large  val¬ 
ues  of  occur  over  the  Gulf  Stream  and  Kuroshio  Systems,  the  Bering  Sea.  the  Arctu. 
Ocean  during  northern  winter  (Fig  5a)  and  over  the  southern  ocean  near  Antarctic  (Fig.  5b) 
since  at  that  time  of  year  prevailing  winds  udvect  cold  air  over  the  waiter  from  neighboring 
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LH  flux  monthly  mean  Jan  1 978 


-30 


-60 
30  60 


Longitude  (E), 


LH  flux  monthly  mean  Jul  1978 


© 

T3 


n 


30 


-60 
30  60 


Longitude 


Fl*t  f>  Monthly  mean  surface  Latent  heat  I'lu*  l  W  mJ  I  simulated  by  theOISS  Coupled  Model  tal 
January  l^'S.  and  tbi  July  I<»1K 

land  masses  ot  ice— sheets. 

Latent  heat  flux 

Figures  6u  and  6b  show  the  GISS  coupled  model  output  of  mean  January  and  July 
(19‘r8)  upward  sensible  heal  flux  at  the  ocean  surface  (0E  ),  respectively.  We  may  find  the 
same  features  as  obtained  from  earlier  studies  such  as  Budyko  (19561:  (a)  The  values  of  Qt 
are  almost  everywhere  greater  in  winter  than  in  summer,  (b)  The  high  values  of  are  found 
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Net  heat  flux  monthly  mean  Jan  1978 


Longitude  (E),  Depth=0(m) 


Net  heat  flux  monthly  mean  Jul  1978 


Longitude  (E), 


-60 
30  60 


Fig  7  Monthly  mean  lurfac*  nel  downward  heal  flu*  IW  m  l  umuUled  hv  the  < ■  ISS  Coupled 
Model  (al  Junuary  Id'H.  and  (b)Jul>  Id78 

near  Ihe  western  boundaries  of  the  North  Atlantic  and  North  Pacific  Oceans  in  northern  win 
ter  (Fig.  6a)  (c)  Seasonal  variations  of  Q E  arc  nowhere  greater  than  over  the  Gulf  Stream 
and  Kuroshio  Systems,  td)  Over  the  eastern  parts  of  oceans  seasonal  variations  are  very  much 
smaller,  because  ut  all  times  of  the  year  prevailing  winds  advect  maritime  air  (moist)  over  the 
cool  ocean  (e)  Throughout  the  year  high  values  of  £)fc  I  TOW  m2)  occur  in  the 
subtropics  and  in  trade— wind  belts, 
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Downward  net  surface  heat  flux 

After  examining  temporal  and  spatial  variations  of  each  component  of  the  surface  energy 
budget,  we  computed  the  net  downward  surface  heat  flux  using  Eq.  (1).  Figures  7a  and  7b 
show  the  GISS  coupled  model  output  of  mean  January  and  July  (1978)  downward  net  heat 
flux  at  the  ocean  surface  ( Q ),  respectively.  It  has  strong  seasonal  variation:  positive  values  in 
summer  hemisphere  and  negative  values  (upward  flux)  in  winter  hemisphere.  Furthermore, 
strong  net  heat  flux  occurs  in  equatorial  regions  with  a  maximum  value  of  140  W  /  m2. 

3.  Statistical  evaluation 


3,1  Cross-correlation  coefficient 


Same  as  in  Chu  et  al.  (1998),  the  validity  of  Haney-type  boundary  conditions  (3),  and  (5) 
should  be  first  tested  by  cross-correlation  coefficients  (CCC)  between  the  two  time  series: 
Q[or  Q]  and  AT  at  day  d,  which  is  computed  in  this  paper  from  a  91 -day  subset  between  45 
days  prior  to  and  45  days  after  that  date, 


Rq, st  (d)~ 


jj  Q{xy)][AT(.x,ysi)-  AJ(x,y)]j 

Oq  (dhAr(d) 


(d-  45  d+  45)  (6) 

where  JV=  91  (day),  (,x,y)  the  horizontal  coordinates,  Q(x,y)  and  A T(x,y)  the  temporal 
means  of  the  subset  data,  and  (d)  and  o  AT  (d)  the  standard  deviations  of  the  subset  data. 

Thus,  we  established  little  over  one  year  (368  days,  17  March  1978—20  March  1979)  CCC 
data  sets:  RQ  AT ,  and  47-  ■  The  higher  the  values  of  CCC,  the  higher  the  confidence  of 
linear  relationship  between  ( Q,Q )  and  AT. 

Besides  the  correlation  analysis,  we  need  to  do  the  variance  analysis,  because  it  does  mat¬ 
ter  if  (Q,Q)  and  AT  correlate  wonderfully,  but  typically  used  values  of  k  end  up  missing  all 
the  variance. 


3.2  Statistical  tests 


Whether  the  sample  values  of  Rs  represent  good  or  bad  linear  relationships  (3),  and 
(5)  between  iQ.Q)  and  AT  should  be  tested.  First,  we  use  the  r— test 


t  = 


V 1  Rqat 


(7) 


to  identify  the  region  with  small  values  of  CCC  in  which  there  is  no  correlation  between  ( Q,Q ) 
and  AT.  We  begin  with  the  null  hypothesis  that  the  two  variables  Q  and  AT  are  not  positively 
correlated.  The  significance  level  a  is  the  probability  that  the  given  value  of  t  is  exceeded  pure¬ 
ly  by  chance.  Equivalently,  it  is  the  probability  of  incorrectly  rejecting  the  null  hypothesis. 
For  a  given  significance  level  a(«=  0.05  was  used),  we  compare  the  r-value  with  ta ,  if  t  <  t, , 
we  accept  the  null  hypothesis,  and  there  is  no  linear  (positive)  relationship  between  Q  and  AT 
at  the  given  location  ( x,y ). 

Second,  we  set  up  a  prior  significant  value  (p0  )  for  CCC,  and  use  the  2-test 
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(1+  RO'AT)(l-p0V 
(1—  ^2.47  Xl+  Po)_ 


(8) 


to  see  if  CCC  is  significantly  larger  than  p0 .  This  z-value  is  satisfying  the  normal  distribution. 

We  start  with  the  null  hypothesis  that  Hq'-R-qat  ~  PoM\ 'Rqat  >  Pa  ■  The  significance 
level  a  is  the  probability  that  the  given  value  of  z  is  exceeded  purely  by  chance.  For  a  given 
significance  level  «(a=  0.05  was  used),  we  compare  the  z-value  with  z„ ,  if  z>  z„ ,  we  accept 
Hl ,  and  the  CCC  between  Q  and  A T  exceeds  p0  at  the  given  location  (x,y).  In  this  study,  we 
usep0  =  0.7. 


4.  CCC  fields  over  the  world  oceans 

The  time  series  RqAT  (x,y,d),  and  RqAT  {x,y,d)  were  decomposed  into  ensemble  means 
^Rqat)  =  Y.Rqat  (atJfitf)  X/Ig.ar)5'  jgg  (x<y>$  >  (9) 

which  show  the  overall  linear  relationships  between  ( Q,Q )  and  AT  and  anomalies 

Rq,st  =  Rq,at~~  >  (10) 

which  indicate  temporal  variations  of  such  relationships.  The  region  with  large  CCC  values 
O  0.7)  is  called  the  valid  region,  and  the  region  with  small  CCC  absolute  values  (<  0,2)  is  cal¬ 
led  the  invalid  region. 

4.1  Ensemble  mean  CCC  fields 

The  ensemble  means  <I?2a7->  (Fig.  8a)  and  Ar  >  (Fig.  8b)  have  a  strong  latitudinal 
variation:  very  small  absolute  values  in  the  equatorial  zone  (10°S-10“N),  and  quite  large  val¬ 
ues  in  the  middle  and  high  latitudes.  These  features  are  very  similar  to  the  analysis  on  the 
NCEP  re-analyzed  surface  data  (Chu  et  aL,  1998).  However,  (Rq  AT  )  and  (Rg^r  )  are  quite 
different  in  the  middle  and  high  latitudes.  The  values  of  (RQAT  >  are  less  than  0.7  almost  ev¬ 
erywhere  except  in  small  regions  of  the  South  Atlantic  Ocean  and  the  western  Indian  Ocean 
between  30°S  and  50°S,  which  is  generally  lower  than  the  analysis  on  the  NCEP  data  (Chu  et 
al.,  1998),  On  the  other  hand,  the  values  of  (Rg  Ar  )  are  greater  than  0.8,  which  is  generally 
higher  than  the  analysis  on  the  NCEP  data  (Chu  et  al.,  1998).  Therefore  according  to  the 
GISS  coupled  model  results,  the  Haney-type  boundary  condition  without  solar  radiation  is  a 
good  representation  in  the  middle  and  high  latitudes  (higher  than  30°S).  The  invalidity  in  the 
equatorial  regions  (10°S-10°N)  may  imply  the  improper  representation  of  latent  heat  flux  QE 
by  temperature  difference  AT.  To  confirm  this,  we  computed  CCC  between  latent  heat  flux 
2e  an<l -  AT  (Fig.  9).  The  world  oceans  are  therefore  divided  into  two  parts:  (a)  equatorial 
zone  with  no  correlation  between  QE  and  —  AT,  and  (b)  raid-  and  high-  latitude  zones  with 
evident  correlation  (>  0.5), 

4.2  Anomaly  CCC  fields 

We  used  Empirical  Orthogonal  Function  (EOF)  analysis  to  delineate  the  major  modes  of 
variability  of  the  anomaly  CCC  field:  R  ’QAT ,  which  is  decomposed  by 

R- 


(11) 
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with  solar  radiation 


longitude  |E).  Depth=0(m) 


wilhout  solai  ladiation 


Longitude  (E).  Doptti  0(m) 


Fi|t  t>  Annual  incun  cTn»»-vomlalir>n  tocfTicicnU  txltt  ten  A/  ami  nr  I  tmlarr  heat  llu\  la  i  with 
solar  nnlulinn  J.  and  ltd  without  solar  radiation  Q 


where  PC  ,td)  is  the  principal  component,  representing  the  temporal  variation  of  the  assoc i 
ated  spatial  pattern  described  by  EOF  p,  ( \  .1  I  The  time  series  is  analogous  to  a  projection 
ofCCC  anomalv  through  the  '  filter  of  an  EOF  mode  during  a  time  scale 

The  first  leading  HOF  is  able  to  account  foroT'o  (7|%)of  the  total  variance  with  (with 
out)  solar  radiation  during  the  integration  period  (Table  I. )  Lach  HOE  mode  is  normalized  so 
that  its  spatial  variance  is  equal  to  units  So.  those  patterns  of  the  first  HOF  arc  enough  to  ex 
plain  the  spatial  anomalies  of  the  global  CCCs.  Here,  we  only  show  EOFJ  and  PC  {  as  an  ex¬ 
ample 
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EOFI  mode  (Fig,  IUa)  has  u  dipole  pattern  featuring  the  Northern  Hemisphere  positive 
(maximum  value  near  0.04)  and  Southern  Hemisphere  negative  (minimum  value  near  — 0  04). 
and  accounts  lor  up  to  67“ »  of  the  total  spatial  variance.  This  pattern  is  believed  to  he  related 
to  the  solar  radiation  The  first  principal  component.  (//>,  during  the  integration  period  is 
shown  in  Fig.  10b.  EOFI  mode  tp ,  (v,  ),  is  generally  positive  (negaiivcl  in  the  Northern 

tSouthern)  Hemisphere  (Fig.  10a).  Therefore.  /,(.'l(t/)>  0.  corresponds  to  positive  (negative) 
CCC  anomalies  in  the  Northern  (Southern!  Hemisphere,  and  PC  (</)<  0,  corresponds  to 
negative  (positive)  CCC  anomalies  in  the  Northern  (Southern!  Hemisphere  From  curb  Oc 
tober  when  PC ,  (t/)~  0,  it  increases  with  time  until  the  end  of  December  when  PC  reaches 
the  maximum  value  of  10.0,  and  then  decreases  to  0  in  late  March.  Alter  late  March.  PC  be 
comes  negative  until  early  October  and  reaches  the  minimum  value  of  -10  0  in  early  July. 
From  early  October  to  early  March,  the  Northern  (Southern)  Hemisphere  is  winter  (summer). 
PC ,  (</)>  0.  The  contribution  of  EOFI  to  the  total  CCC  field.  PC  { it/)*  «»,  (.v  ,v, ),  adds  pos 
Hive  values  to  the  ensemble  CCCs  with  a  maximum  increasing  value  of  0.41  ~  10.0*  0,041  in 
the  Northern  Hemisphere  and  adds  negative  values  to  the  ensemble  CCCs  with  a  maximum 
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Fig.  10  <  a  I  F  ii  »t  EOF  mode  of  H  .  Al  .and  Hit  *i»ocmied  tuiK  «iif*  of'/’C, 

decreasing  value  of  -0.4[  =  10.0*  (—0,04)1  in  the  Southern  Hemisphere.  On  the  other  hand, 
from  early  March  to  early  October.  PC  { (</)<  0  The  contribution  of  FOF1  to  the  total  C'CC 
field  is  reversed  to  the  period  from  early  October  to  curly  March.  Thus,  the  winter  hemisphere 
has  better  cross-correlations  than  the  summer  hemisphere  The  maximum  temporal  variabili 
tv  can  reach  0,8  (  =  2  *  0.4>  Thus.  Ihe  Haney  type  thermal  boundary  condition  with  solar  ra¬ 
diation  is  valid  only  in  the  winter. 

The  physical  interpretation  to  Ihe  results  that  the  winter  hemisphere  has  better 
cross-correlation  between  heat  (lux  and  air-sea  temperature  difference  than  the  summer 
hemisphere  is  that  the  values  of  /<s  are  almost  everywhere  smaller  in  winter  than  in  summer, 
and  /?  s  cannot  be  parameterized  by  the  air-sea  temperature  difference 

5.  Seasonal  C'CC  variation  in  five  geographic  /ones 

On  the  busis  of  the  spatial  variation  of  the  ensemble  mean  field,  we  divided  Ihe  world 
oceans  into  live  geographic  zones;  equatorial  zone  ( HflM-KTS).  northern  subtropical  zone 
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Fig.  1 1.  Seasonal  variation  of  cross— correlation  coefficients  averaged  within  five  geographic  zones 
between  (a)  Q  and  AT  ,  and  (b)  and  AT.  Here,  the  thin  solid  curve  indicates  the  southern  middle 
and  high  latitude  zone  (south  to  30°S),  the  thick  solid  curve  indicates  the  northern  middle  and  high 
latitude  zone  (north  to  30°N),  the  dotted  curve  indicates  the  southern  subtropical  zone  (10-30°S), 
the  dashed-dotted  curve  indicates  the  northern  subtropical  zone  (10-30°N),  and  the  dashed  curve 
indicates  the  equatorial  zone  (10°S— 10°N). 
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(10-30°N),  southern  subtropical  zone  (10— 30“S),  northern  middle  and  high  latitude  zone 
(north  to  30°N),  and  southern  middle  and  high  latitude  zone  (south  to  30°S).  At  each  day,  we 
averaged  the  CCC  data  RQ  iT  (x,y,d),  and  R^  at  (x,y,d)  spatially  within  each  zone.  Thus,  we 
obtained  five  time  series  for  each  of  the  CCC  data;  R^  a7-(d),^  Ar  (d).  Here,  the  superscript 
‘  S'  means  the  spatial  average. 

5.1  W  ith  solar  radiation 

The  seasonal  variation  of  AT id)  for  the  five  geographic  zones  (Fig.  11a)  is  much 
stronger  than  that  for  the  NCEP  data  (Chu  et  al.,  1998).  The  value  of  Rq  at  (d)  fluctuates 
from  -0.25  to  -0.08  in  the  equatorial  zbne  (no  correlation  between  Q  and  A7V,  from  0.04  to 
0.54  in  the  southern  subtropical  zone,  and  from  -0.2  to  0.5  in  the  northern  subtropical  zone 
(weak  correlation  between  Q  and  AT).  However,  the  cross— correlation  coefficient  Rq  at  W) 
has  a  seasonal  variation  (0.2  to  0.86)  in  both  the  northern  and  southern  middle  and  high  lati¬ 
tude  zones  with  180°  phase  shift:  higher  values  during  winter  and  lower  values 
during  summer.  Such  a  feature  (winter  hemisphere  has  a  higher  correlation)  is  due  to  small 
values  of  solar  radiation  (Rs )  in  winter  hemisphere.  Among  the  four  components  in  the  net 
surface  heat  flux  (Q),  only  solar  radiation  cannot  be  represented  the  air-sea  temperature  dif¬ 
ference  AT. 

The  period  when  R |  (d)>  0.7  (valid  period)  is  from  25  September  to  15  March  in  the 

northern  middle  and  high  latitudes  and  from  30  March  to  1  September  in  the  southern  middle 
and  high  latitudes.  Such  a  strong  seasonal  variation  of  R^  AT  (d)  in  the  middle  and  high  lati¬ 
tudes  indicates  the  improper  representation  of  the  solar  radiation  Rs  by  the  temperature  dif¬ 
ference  AT. 

5.2  Without  solar  radiation 

Figure  lib  illustrates  the  seasonal  variation  of  8gi47.((()  for  the  five  geographic  zones. 
The  correlation  is  higher  between  Q  and  AT  than  between  Q  and  AT,  indicating  the  benefit 
of  excluding  the  solar  radiation  from  the  net  surface  heat  flux.  For  example,  id)  is  near 
0.86  all  year  round  in  the  southern  middle  and  high  latitudes,  and  varies  between  0.74  and 
0.88  in  the  northern  middle  and  high  latitudes.  However,  R^At  (d)  is  still  very  low  (between 
-0,06  and  0.06)  in  the  equatorial  zone,  and  becomes  higher  but  less  than  0.7  in  both  northern 
and  southern  subtropical  zones. 


6.  Seasonal  variation  of  relaxation  coefficient  k 


Since  the  value  of  k  is  quite  uncertain,  such  a  variance  comparison  between  (Q  ,Q )  and 
AT  depends  on  the  choice  of  k.  If  k  is  determined  by 


x(d)= 


,(d) 


-(d) 


'g.sr 


(12) 


at  each  time  instance  d,  the  variance  of  Q  equals  the  variance  of  kAT.  This  is  so-called 
no— missing  variance  to  estimate  Q  by  kAT.  Now,  we  should  check  if  values  of  k  calculated 
by  (12)  are  within  the  reasonable  range. 
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Fig.  12.  Seasonal  variation  of  relaxation  coefficient  k  averaged  within  five  geographic  zones,  de 
termined  by  no  loss  of  variance  for  (a)  Q  =  kAT,  and  (b)  Q  =  kAT,  Here,  the  thin  solid  curve  indi¬ 
cates  the  southern  middle  and  high  latitude  zone  (south  to  30°N),  the  thick  solid  curve  indicates  the 
northern  middle  and  high  latitude  zone  (north  to  30°S),  the  dotted  curve  indicates  the  southern 
subtropical  zone  (10-30°S),  the  dashed-dotted  curve  indicates  the  northern  subtropical  zone 
(10-30°N),  and  the  dashed  curve  indicates  the  equatorial  zone  (10°S-10°N). 
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We  computed  the  relaxation  coefficient  k  using  (12)  and  then  averaged  k  spatially  within 
each  of  the  five  zones.  Thus,  we  obtained  both  the  spatial  (five  zones)  and  the  temporal  varia¬ 
bility  of  k.  Theoretically,  k  can  be  time-dependent.  Practically,  k  is  usually  taken  as  a  con¬ 
stant  in  the  surface  thermal  boundary  conditions  (3),  and  (5),  Therefore  besides  high  cross 
correlations,  we  should  also  consider  quasi-steadiness  of  k  in  evaluating  the  surface  thermal 
boundary  conditions. 

For  surface  net  heat  flux  with  solar  radiation  (Fig.  12a),  k  varies  from  —90  to  —20 
Wm_2K_1  in  the  equatorial  zone,  from  0  to  130  Wm't1  in  the  southern  subtropical  zone, 
from  -20  to  124  Wm"2K_1  in  the  northern  subtropical  zone,  from  20  to  100  Wm~!K_1  in  the 
southern  middle  and  high  latitude  zone,  and  from  20  to  120  Wm"2K"'  in  the  northern  middle 
and  high  latitude  zone.  When  the  solar  radiation  is  excluded  from  the  net  surface  heat  flux 
(Fig.  12b),  k  has  small  temporal  variations  in  all  the  five  geographic  zones.  For  example,  k 
varies  from  80  to  120  Wm'f  in  the  northern  middle  and  high  latitudes  and  front  100  to  1 20 
Wm_2K_1  in  the  southern  middle  and  high  latitudes. 

On  the  basis  of  (a)  high  correlation  between  Q  and  A Tx ,  and  (b)  quasi-steadiness  of  k  , 
we  may  conclude  that  the  surface  thermal  boundary  condition 

5  =  kAT 

is  a  good  parameterization  of  surface  net  heat  flux  for  the  middle  and  high  latitudes. 

7.  Variance  analysis  for  the  haney-type  condition 

The  coupling  coefficient  k  has  units  of  watts  per  square  meter  per  kelvin.  A  time-scale 
can  be  derived  as  r=  pcp  h  /  k,  where  p  is  the  density,  cp  is  the  specific  heat  of  seawater,  and 
h  is  a  typical  mixed  layer  depth.  Weaver  and  Sarachik  (1991)  used  a  time— scale  of  25  days. 
Marotzke  and  Stone  (1995)  used  30  days.  For  a  mixed  layer  50  m  deep  these  correspond  to 
values  of  k  of  97  Wm_2K_1  and  81  Wm't1,  respectively.  As  pointed  by  a  number  or  au¬ 
thors  (Rahmstorf  and  Willibrand,  1995;  Pierce  et  al,,  1995;  Cai  and  Chu  1996),  change  of  val¬ 
ues  of  k  has  detrimental  effects  on  the  realism  of  a  modeled  thermohaline  circulation.  Thus, 
use  of  a  proper  k- value  becomes  an  important  issue. 

We  computed  the  standard  deviations  of  Q,  and  kAT  averaged  over  the  northern  and 
southern  middle  and  high  latitude  zones:  aQ  (d)  and  <rKi7-  (d).  The  relaxation  coefficient  k  has 
values  of  30,  60,  90,  and  120  Wm'V.  For  a  given  day,  d,  the  standard  deviation  aKAT  (d)  in¬ 
creases  with  k.  When  k  =  120  Wm‘!K"',  the  two  standard  deviations,  og  (d\  and  o (d),  are 
quite  close  in  both  northern  (Fig.  13a)  and  southern  (Fig.  13b)  middle  and  high  latitude 
zones.  The  large  value  of  k  (120  Wm'f)  might  be  caused  by  a  short  coupling  time-scale, 
since  in  the  GISS  coupled  model  atmospheric  and  oceanic  surface  fluxes  of  water,  heat  (ex¬ 
cluding  solar  radiation),  and  momentum  are  of  opposite  sign  and  are  applied  synchronously. 

8.  Conclusions 

On  the  basis  of  cross-correlation  and  variance  analyses,  we  obtain  the  following  results: 

(1)  The  Haney-type  conditions,  Q=  kAT  and  Q  =  kAT,  have  better  representation  of 
surface  heat  forcing  in  the  middle  and  high  latitudes  than  in  the  equatorial  regions.  On  the  ba¬ 
sis  of  high  correlation  between  Q  and  AT  and  associating  quasi-steadiness  of  k. 
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Fig,  13.  Comparison  among  ag  (dashed  curve),  and  at  different  values  of  tc(A:s:=  30 
Wm'!KH,  B:  k=60  Wm"2K"',  C:  k=»0  Wm'3K  ',  and  D:  »c=  120  WnT2K"')  for  (a)  northern 
middle  and  high  latitudes,  and  (b)  southern  middle  and  high  latitudes. 
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we  recommend  to  use  0=  kAT  for  the  middle  and  high  latitude  surface  thermal  forcing. 

(2)  For  the  equatorial  and  subtropical  oceans,  the  Haney-type  conditions,  Q  =  kAT 
and  Q  =  kAT,  are  not  good  approximations  for  the  surface  thermal  forcing. 

(3)  A  value  of  120  Wm-2K_1  for  the  relaxation  coefficient  k  is  suggested.  When  k  is 
taken  as  120  Wra"!K_1,  the  variance  of  surface  net  heat  flux  Q  can  be  well  represented  by  the 
variance  of  kAT.  This  value  is  much  higher  than  that  estimated  from  the  NCEP  data. 

(4)  Quantitative  values  of  correlation  coefficient  and  relaxation  coefficient  may  vary  with 
the  eddy  viscosity  and  diffusivity.  But  the  results  on  validity  of  the  Haney— type  boundary 
condition  holds. 

Authors  are  indebted  to  Gary  L.  Russell  at  NASA/  Goddard  Institute  for  Space  Studies  for  allowing  us  to  use 
the  coupled  atmosphere-ocean  model.  This  work  was  funded  by  the  Office  of  Naval  Research,  Naval  Oceanographic 
Office,  and  Naval  Postgraduate  School 
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